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A theoretical model for studying pattern formation in electroconvection is proposed in the form 
of a modified Swift-Hohenberg equation. A localized state is found in two dimension, in agreement 
with the experimentally observed “worm” state. The corresponding one dimensional model is also 
studied, and a novel stationary localized state due to nonadiabatic effect is found. The existence of 
the ID localized state is shown to be responsible for the formation of the two dimensional “worm” 
state in our model. 
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When a system is driven away from its equilibrium 
state, it often responds in forming regular patterns. Pat¬ 
tern formation phenomena is rather ubiquitous, occur¬ 
ring in many different physical, chemical and biological 
systems. One of the most carefully studied systems is the 
Rayleigh-Benard convection (RBC), where high precision 
experimental measurements can be carried out and com¬ 
pared quantitatively with theoretical studies. The care¬ 
ful study of this rich physical system for the past twenty 
years has greatly advanced our understanding of systems 
far away from equilibrium ]j]]. 

The study of localized structure in nonequilibrium sys¬ 
tem has received a great deal of attention since being ob¬ 
served experimentally in binary-mixture RBC |2j. Even 
though binary-mixture RBC is a highly dissipative sys¬ 
tem, the localized structures behave much like soliton in 
integrable systems. For example, in certain parameter 
range, they can pass through each other without chang¬ 
ing their structures. On the theory side, Thual and Fauve 
H were the first to study the behavior of a subcritical 
complex Ginzburg-Landau equation and found that in 
certain parameter range, there are indeed localized pulse 
solutions. The basic ingredients for the existence of lo¬ 
calized structure are: l)there has to be linear bistability, 
which guarantees the local stability of the peak and the 
tail of the pulse; 2)the nonlinear dispersion (the complex 
part of the coefficients for the nonlinear terms) is needed 
to stabilize the front connecting the peak and the tail of 
the localized solution. Both of these requirements seem 
to be consistent with the experimental system, because 
the initial bifurcation in binary-mixture RBC is subcriti¬ 
cal Hopf bifurcation. Much work has since been devoted 
along these lines to understand the details of the experi¬ 
mental results (|||. 

Most of the experimental results in binary-mixture 
RBC were obtained in quasi one dimension, i.e., in a 
thin annulus (b|. Further efforts to extend these findings 
to two dimensional system have not revealed any similar 
2D localized state as in ID except for some time depen¬ 
dent patchy structure |7| and some long time transients 
H. Recently, M. Denning et al || studied electroconvec¬ 
tion in nematic liquid crystal very carefully. Due to the 
anisotropic nature of the liquid crystal, they found that 


the initial instability of the system is towards forming 
oblique rolls with certain angle with respect to the direc¬ 
tor of the nematic liquid crystal. They also found that 
the initial bifurcation is Hopf bifurcation jnj. Depending 
on the electrical conductivity, the pattern they observed 
above onset is either spatially extended time dependent 
state/spatial-temporal chaos state(STC) or some isolated 
localized state, which they named the “worm” state. The 
worm state is localized in the direction perpendicular to 
the director of the liquid crystal, but is extensive in the 
parallel direction. The worm state can move in the par¬ 
allel direction. The internal structure of the worm state 
seems to consist of both orientations of the oblique rolls, 
and internal roll structure is moving relative to the mo¬ 
tion of its envelop. It is the goal of this paper to under¬ 
stand the interesting structure of the worm state. 

There are a few well established methods for studying 
pattern formation. Notably among them is the ampli¬ 
tude equation formalism which was used in [3] and related 
works HH] to study the pulse pattern in binary-mixture 
RBC. Amplitude equation describes the large scale and 
long time behavior of the envelop of the pattern. It is per¬ 
turbative in nature and describes the system accurately 
at parameters close to the onset. The shortcoming of 
the amplitude equation is thus quite obvious. In writing 
down the amplitude equation, one has already broken the 
full spatial symmetry of the original system which might 
be important in 2D. In the case of subcritical bifurcation, 
the amplitude equation can only give qualitative results 
because the amplitude does not scale with the small pa- 
rameter(reduced Rayleigh number). For the problem at 
hand, there are even more severe limitation for the am¬ 
plitude equation. Since the spatial extension of the worm 
state in the perpendicular direction is only about a cou¬ 
ple of wavelengths, there is no seperation of length scales 
to justify the use of the amplitude equation. 

To understand the experiments qualitatively, a phe¬ 
nomenological model is often useful. The Swift- 
Hohenberg(SH) equation |ll) is an order parameter equa¬ 
tion with the full symmetry of the original problem, and 
its linear properties agree with that of the original prob¬ 
lem. The SH equation is phenomenological in nature, 
and usually can not be derived from the original equa- 
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tions. For different experimental systems, The Swift- 
Hohenberg equation can be modified in different ways 
from its original form to accommodate different physi¬ 
cal situations and symmetry requirements, e.g., the Non- 
Bousinesq effect, mean flow effect 0 ,Hopf bifurcation 
00 and rotating convection Jl5|]. Because of its versa¬ 
tility and simplicity, the modified Swift-Hohenberg equa- 
tion(MSHE) has become instrumental in understanding 
many pattern forming systems. 

In writing down the modified Swift-Hohenberg equa¬ 
tion for electroconvection, we know that the equation has 
to be anisotropic even at the linear level, and the equa¬ 
tion has to be complex because the initial bifurcation is 
Hopf bifurcation. Let <p(x, t ) be the complex order pa¬ 
rameter, we can write the order parameter equation as: 

d(j)/dt = (e + iu)(p - <r((d 2 + q 2 ) 2 + b(d 2 + q 2 x )(d 2 y + q%) 

+ ( d 2 + q 2 y ) 2 )(t> + iv g ((d 2 + q 2 x ) + a(d 2 + q 2 y ))cp 

+ 9o\<P\ 2 (p + gi\<P\ 4 <P (1) 

e is the reduced Rayleigh number, w is the Hopf fre¬ 
quency, q = (■ q x , q y ) = ( cos9 , sinO ) is the linearly most 
unstable wavevector, b is an anisotropic parameter with 
the constrain |6| < 2, and a is a complex constant. The 
first two lines on the RHS of eq. (1) represent the lin¬ 
ear properties of the elctroconvection system. It is easy 
to see that the system is most unstable at \k x \ = q x 
and \k y \ = q y for cp ~ exp(ik x x + ik v y). v g is propor¬ 
tional to the group velocity and a is another anisotropic 
parameter. When a = 1, the group velocity is along 
the wavevector direction q. The last line on the RHS 
of eq. (1) contains the nonlinear coupling terms with 
complex coefficients go and g\. In principle, the nonlin¬ 
ear terms can also be anisotropic, we only include the 
simplest terms possible here. 

Since it is the goal of this paper to find the localized 
worm state, we focus our attention on the subcritical case 
where e < 0, Re(go) > 0 and Re(g\) < 0. We can easily 
get rid of the iui term in the linear part of the equa¬ 
tion by a change of variable cp = e lut (p, so we will set 
ui = 0 for now on. There are five real parameters: e, a, 
b , 9 and v g and three complex parameters: cr, go and g\ 
for this model. We have numerically studied the MSHE 
extensively in the parameter space and identified certain 
parameter region where localized worm state is observed. 

To demonstrate the existence of the worm state, we 
first show the behavior of the eq. (1) for a particular set 
of parameters: e = —0.2, a = 1, b = 0, 9 = 23°, v g = 0.5, 
cr = 1.5, go = 3 + i and g\ = —2.75 + i. The equa¬ 
tion is simulated in systems of size 64 x 64, 128 x 64 and 
256 x 64 with periodic boundary condition using both 
the second order finite difference method and spectral 
method with discretization Aa; = Ay = 0.5,1.0 and time 
step At = 0.001,0.01. We start the system with ran¬ 
dom initial condition with large enough amplitude. The 
system quickly organizes itself into the worm like state. 
A snap shot of the 2D pattern for Re{<p{x,y)) after the 
initial transient is shown in figure la. 


A slice of the 2D pattern along the y direction at x = 45 
is shown in figure lb. The localization of the worm states 
in the y direction is apparent from fig. lb. The Worm 
states travel in the x direction. From their length, the 
worm states in our simulation can be divided into two 
categories, which we call long worm and short worm. The 
length of the short worm does not change with time, and 
is usually ~ 20, which is 3 basic wavelengths long. They 
travel in the x direction with constant velocity propor¬ 
tional to v g . An example of a short worm can be seen 
near the bottom of fig.la. The long worm’s length grows 
with time and eventually extends over the whole length of 
the system because of the periodic boundary condition. 

We have tested the sensitivity of the worm pattern 
to the parameters in our model. We find that there is 
a finite range of parameters where the worms appear. 
For example, if we change the value of e while keeping 
the rest of the parameter unchanged, worm exists for 
—0.10 > e > —0.25. When e is too small, there is no 
pattern; and when e is too big, the pattern becomes ex¬ 
tended instead. The worm state is quite insensitive to 
the values of a and b , as long as a ~ 1 and |6| < 2. 
For 6 = 2 and a = 1, the model becomes isotropic and 
the worm structure gives away to time dependent patchy 
structure 0]. The velocity v g is important to give the 
worm a group velocity. The wavevector angle 9 has to 
be small enough 9 < 35° to make the worm perfectly 
aligned in the x direction. There are also finite regions in 
the parameters go , <?i, cr where worm states are observed. 

The worm states interact strongly with each other. 
When two short worms collide, they come out of the 
collision without changing their characteristics. When 
a short worm collides with a long worm, the short worm 
sometimes disappears. When two long worms approach 
each other off center, oblique rolls are excited in the re¬ 
gion of their overlap until the worms pass through each 
other or one of the worms disappears. When two long 
worms collide head on, they stop each other and form a 
well defined boundary between them. 

For the short worm, because of the spatial extention in 
both directions are about the same order, the formation 
of the short worm should be due to strong interaction be¬ 
tween the two dimensions. However, for the long worms, 
due to the extensiveness of the worm in the x direction, 
we are able to separate the dependence in the two di¬ 
mensions and therefore gain more understanding of the 
mechanism for the localization in the y direction. In¬ 
deed, Fourier analysis of the long worm along x direction 
shows that it is a good approximation to assume the x 
dependence to be a simple plane wave: 

<P(x, y, t ) = ip(y, t) exp(ik x x) (2) 

If we substitute the above ansatz into the original equa¬ 
tion (1), we obtain a one dimensional dynamical equation 
for ip(y, t). For simplicity, we set a = 1 and 6 = 0: 

dip/dt = (e + iu)ip - a{d 2 y + q 2 y ) 2 ip + iv g {d 2 + q 2 y )ip 

+ golV’l 2 ^+ 5i|V’lV (3) 
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where, e = e — Re(a)(q% — k ^) 2 and Co = to — Im(a)(ql — 
k^.) 2 + v g(Qx — k 2 )- For the same reason as in two dimen¬ 
sion, we can set Co = 0. 

We have studied the above ID MSHE carefully. The 
numerical scheme is the same as in the two dimensional 
case, and we also started with a random initial condi¬ 
tion with sufficient amplitude. In order to compare it 
to the two dimensional case, we have set the parameters 
a = 1.5, v g = 0.5, go = 3.0 + i and g\ = —2.75 + i to be 
the same as in the 2d calculation. We can vary e because 
the value of k x is undetermined a priori. For q y = sin(6) 
with 9 = 23°, we found a finite range of e values, where 
localized structure is observed —0.15 > e > —0.5. A 
space-time plot of Re(i/>(y,t )) is shown in figure 2. It is 
clear that the system evolves to a final state with two 
localized pulses. Most remarkably, the pulses are not 
moving even in the presence of the group velocity term 
in eq. (3). 

We found that the pulse solution can be written as: 

i>(y, t ) = A{y) exp (ia{y, t)) (4) 

where the amplitude A(y) is independent of time and is 
localized with a width of ~ 10. Shifting the peak posi¬ 
tion to be at y = 0, the shape of the pulse is symmetric 
around y=0: A{y) = A{—y) . The phase of the pulse 
depends on time linearly: 

a{y,t) = a 0 (y)+ M (5) 

with H = —.25. The shape of the time independent 
phase ao(y) is depicted in figure 3(b). From figure 
3(b), we see that the phase is symmetric around y = 0: 
cno(y) = ao(—y). Away from the peak of the pulse, the 
phase is roughly linear in y: 

oto{y) ~ k v \y\ + const \y\ > 5 

with k y ~ 1. 

As we pointed out earlier in our paper, the existence 
of localized state in subcritical equation with complex 
coefficient was known |j). However, a stationary local¬ 
ized pulse in the full equation including group velocity 
term is observed for the first time here. If one were able 
to eliminate the small scale structure and write the full 
equation in terms of the amplitude equation, one could 
use two amplitudes characterizing the left and the right 
moving wave packets. As shown in the work of Brand 
and Diessler [^6|, the left and the right moving pulses 
often can pass through each other without altering their 
own characteristics. Upon tuning the inter-coupling be¬ 
tween the left and the right moving pulses, the pulses 
can form bound state which do not move in either direc¬ 
tions. However, the structure of the bound state is such 
that the amplitudes of both the left and the right moving 
pulses are strongly suppressed in their coexisting region. 
In the middle of the bound state, both amplitudes are 
very small, and overall, within the bound state, both 
pulses still keep their own identities. 


The pulse structure we observed here can not be ex¬ 
plained as the bound state of the left and the right mov¬ 
ing pulses in the amplitude equation because the local¬ 
ized state here does not have local minimum at the cen¬ 
ter. On the contrary, the amplitude is the maximum at 
the center of the localized pulse. Furthermore, in our sim¬ 
ulation, we never observe any individual traveling pulse 
and the stationary localized structure always forms spon¬ 
taneously as an whole object. The two halves of the pulse 
seem to have opposite phase velocity v p = ±Q/k y . How¬ 
ever, the periodicity in the phase 2 Tt/k y is much smaller 
than the linearly most unstable wavelength 2ir/q y . In 
fact, the size of the pulse, i.e., the spatial extension of 
the whole pulse is smaller than 2n/q y . This clearly shows 
that the amplitude equation approach is invalid here. 
Therefore we believe that the localized pulse observed 
here is indeed a novel structure due to non-adiabatic ef¬ 
fect, which can only be studied using MSHE where vari¬ 
ation in all the length scales are kept. 

According to the above analysis, the long worm can be 
understood as the combined structure of the Id localized 
state in y direction and the simple harmonic behavior 
in x direction. In the y direction, the amplitude of the 
worm is maximum at the center and decays symmetri¬ 
cally as we move away from the center. The phases on 
the two sides have effective wavevectors k± = (k x ,±k y ) 
which are different from the linearly most unstable modes 
q± = (q x ,±q y ). The long worm expands while moving 
with the group velocity. 

In summary, we have proposed a modified Swift Ho- 
henberg model to count for the formation of localized 
worm state seen experimentally in electroconvection ex¬ 
periment. For a broad parameter range, we have found a 
solution of the MSHE which is localized in one direction 
and extended in the other. The structure of our solution 
resembles the experimental findings closely. The local¬ 
ization of the solution is further understood by simplify¬ 
ing the original two dimensional equation to ID. A novel 
localized stationary pulse state is discovered in the ID 
study which explains the worm structure. 

To understand the experiments fully, we need to un¬ 
derstand the dynamics of the worm in detail. For our 
current model, once a worm is formed, it will persist until 
it interacts with the boundary or another worm, while in 
the experiment, it seems that it can dissappear by itself. 
Such a phenomena could be due to the strong external 
fluctuation or certain internal instability, which depends 
on the details of nonlinear terms. Another interesting is¬ 
sue is the interaction between worms, which also depends 
crucially on the form of the nonlinear terms. The MSHE 
can also be used to study the spatially extended STC 
state observed for larger electrical conductivity and the 
transition between the worm state and the STC state. 

The author would like to acknowledge useful discussion 
with Drs. M. C. Cross , H. Riecke and M. Dennin. 
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FIG. 1. (a)Grey scale plot of Re{<j>{x,y)) showing localized 
worm structures. (b)The profile of the 2D pattern in (a) at 
x=45: Re((j)(4:5, y)) versus y, the localized nature of the worm 
state is obvious. 

FIG. 2. Space time plot of Re(ip(x,t)) + i/4 versus x for 
time difference dt = 4 

FIG. 3. (a)The amplitude of the ID pulse A(y) versus y\ 
(b)The stationary part of the phase of the ID pulse ao(y) 
versus y. 
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